Tunable Surface Conductivity in Bi 2 Se 3 Revealed in Diffusive Electron Transport 
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We demonstrate that the weak antilocalization effect can serve as a convenient method for detect- 
ing decoupled surface transport in topological insulator thin films. In the regime where a bulk Fermi 
surface coexists with the surface states, the low field magnetoconductivity is described well by the 
Hikami-Larkin-Nagaoka equation for single component transport of non-interacting electrons. When 
the electron density is lowered, the magnetotransport behavior deviates from the single component 
description and strong evidence is found for independent conducting channels at the bottom and 
top surfaces. Magnetic-field-dependent part of corrections to conductivity due to electron-electron 
interactions is shown to be negligible for the fields relevant to weak antilocalization. 



The surface of a 3D topological insulator (TI) [l], Q 
hosts a 2D system of Dirac electrons with spins trans- 
versely locked to their translational momenta. Such spin- 
helical surface states Q offer a new route for realizing 
exotic entities such as Majorana fermions and magnetic 
monopoles Q ■ The unique surface spin structure also has 
profound impact on the transport properties of TI [H [|| . 
The Berry phase associated with the surface electrons 
causes suppression of backscattering Q and hence im- 
munity to localization regardless of the strength of disor- 
der [5j. This weak antilocalization effect can be brought 
out by applying a perpendicular magnetic field. It pro- 
duces a negative magnetoconductivity due to the break- 
ing of time-reversal symmetry. The negative magneto- 
conductivity has indeed been observed in various TI thin 
films by several groups [7-ll|. However, most of these 
measurements were carried out with samples in which 
the Fermi energy is not located in the band gap, so 
that they are not in the so-called topological transport 
regime [lj; 



13(| . Since topologically trivial 2D electron 



systems (e.g. Au thin films) may also exhibit similar mag- 
netoconductivity behavior as long as the spin-orbit cou- 
pling (SOC) is sufficiently strong concern has to be 
raised whether weak antilocalization can provide a reli- 
able method for identifying the surface state transport, 
which is a key starting point for future exploration of 
various topological effects and novel devices [1, 15- 19j. 

Here we confirm unequivocally that the weak antilo- 
calization effect can be used to differentiate the surface 
transport from transport dominated by bulk carriers. 
This is demonstrated on Bi2Se3 thin films with carrier 
densities that can be tuned over a wide range with a 
back-gate. When the transport is not in the topolog- 
ical regime, the magnetoconductivity can be described 
by a single component Hikami-Larkin-Nagaoka (HLN) 
equation (20[. This description is found to be valid 
for a remarkably wide range of electron densities (0.8- 
8.6xl0 13 cm~ 2 ) in samples with the Fermi energy located 
inside the conduction band, even if electron-electron in- 
teractions are taken into account. In contrast, in a regime 



where the electronic system is split up into an electron 
layer at the top surface and a hole layer at the bot- 
tom, the magnetoconductivity deviates strongly from the 
single-component HLN equation. Our analysis provides a 
convenient method for detecting decoupled surface trans- 
port. It complements existing techniques based on quan- 
tum oscillations that are limited to samples of high car- 
rier mobilities 
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or samples with a quasi-lD geome- 
try [22j. 

The Bi2Se3 thin films were grown on SrTiC>3(lll) sub- 
strates with molecular beam epitaxy [23j]. The dielectric 
properties of SrTi03 are well suited for gating purposes 
and the carrier density in these devices can be varied by 
at least 2 x 10 13 cm~ 2 [7|. All of the samples used in 
this work were patterned into 50 /im wide Hall bars with 
photolithography, followed by Ar plasma etching (Fig. 1 
inset). This eliminates uncertainties in evaluating resis- 
tivities encountered in previous transport studies due to 
the influence of electrical contacts or the irregular shape 
of the sample. A set of more than ten samples with thick- 
nesses between 5 and 20 nm has been measured. Most of 
the samples have a back-gate deposited at the bottom of 
the substrate, and a few of them are equipped in addition 
with a top-gate. The latter was deposited on an AlO^ 
layer prepared with atomic layer deposition. Transport 
measurements were carried out in cryostats with temper- 
atures as low as 10 mK and magnetic fields up to 18 T. 

Fig. 1 displays typical magnetotransport data. All of 
the samples show a positive magneto-resistance with a 
sharp cusp around zero magnetic field, consistent with 
previous measurements of Bi2Se3 thin films [7Hllj. As 
demonstrated in Fig. 1(b), the low field magnetoconduc- 
tivity, defined as A<r(B) — u xyi {B) — <7xx(0), can be fitted 
well with the HLN equation in the strong SOC limit, i.e. 
when > r so , r e : 



Aa(B) ~ -a ■ — 
Tin 



(1) 



where t so (r e ) is the spin-orbit (elastic) scattering time, 
if) is the digamma function, B^ — h/(4DeT^,) is a char- 
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FIG. 1. (color online) (a) Magnetoresistance (MR, defined as 
Pxx(B) / Pxx(0) — 1) and Hall resistance of a typical Bi2Se3 thin 
film at T = 1.2 K. The low field MR (between ±3 T) is shown 
more clearly in the lower inset. The upper inset is an optical 
image of the Hall bar used for the transport measurements, 
(b) Magnetoconductivity data (symbols) from four samples 
fitted with Eq. (1) (lines). The dephasing field B$ varies by 
nearly a factor of 30. Extracted values of a for 10 samples 
are plotted as a function of electron density n e and mobility 
(j, in (c) and (d), respectively. Both are evaluated based on 
the low field transport measurements. 



acteristic field related to the dephasing time r^, D is the 
diffusion constant and h is the Planck constant. The coef- 
ficient a takes a value of 1/2 for a traditional 2D electron 
system with strong spin-orbit coupling. The same value 
is expected for the electron transport on one surface of a 
3D TI with a single Dirac-cone (jf . 

Fig. 1(c) shows that the extracted a values are dis- 
tributed in a narrow range near 1/2 for 10 samples 
with 2D electron densities n e spreading from 0.8 to 
8.6xl0 13 cm~ 2 [24j |. No correlation is found between a 
and the electron mobility fj,, which varies nearly two or- 
ders of magnitude (Fig. 1(d)). Based on angle-resolved 
photoemission measurements [ijj], the top and bottom 
surfaces of a Bi2Se3 thin film can only accommodate a 
total electron density of ~ 0.5 x 10 13 cm~ 2 even if the 
Fermi energy reaches the bottom of the conduction band. 
Thus we anticipate a significant number of bulk electrons 
(or quasi-2D electrons with parabolic dispersion) for the 
above range of n e . The nonlinear Hall resistivity curves 
(see e.g. Fig. 1(a)) also suggest the coexistence of multi- 
ple charge carrier types. Even if so, the analysis of the 
weak antilocalization effect itself at small magnetic fields 
yields values of a close to 1/2. In this magnetic field 
regime where the antilocalization effect is observed, these 
samples do behave like 2D systems with a single type of 
charge carrier. This can only be understood when there is 
a strong mixing between the surface and the bulk electron 
states or when the dephasing field of one of the conduct- 




FIG. 2. (color online) (a) Magnetoconductivity Aa(B) at T = 
2 K and Vg = plotted as a function of perpendicular magnetic 
field B± for several tilt angles. 6 — refers to B perpendicular 
to the thin film plane, (b) Temperature dependencies of (T X x 
(open symbols) recorded at B =0, 0.2, 1, and 5 T. The straight 
lines are linear fits of ct X x to InT. In the upper inset, the 
slope, defined as k — (irh/e 2 )dax X (B,T)/d(lnT), is plotted 
as a function B. The electron density is about 2 x 10 13 cm -2 
so that Ef is located in the conduction band. 



ing components (i.e. the bulk or top/bottom surface) is 
much smaller than those of the others. We note that it 
was demonstrated long ago that the two-valley 2D elec- 
tron system confined in a Si inversion layer displays a val- 
ues close to that for a single- valley system and not the one 
expected for two independent valleys |25j . Fukuyama at- 
tributed it to intervalley scattering [26j. Similar physics 
might take place here because of considerable scattering 
between the surface and bulk states when the Fermi en- 
ergy is located in the conduction band. 

The robustness of a ~ 1/2 is at first sight surpris- 
ing, because sources other than weak antilocalization 
may also contribute to the low field magnetoconductiv- 
ity. In the non-interacting electron picture, the Zeeman 
energy, which was not considered in deriving the HLN 
equation, is known for mixing the spin singlet and triplet 
states and hence suppressing the weak antilocalization 
effect |27j. The corresponding correction to Aa-(B) is 
determined by the ratio 7 = Ez/E so — g^isB / (hr^ 1 ), 
where g is the electron g-factor. The Zeeman energy 
also causes an extra change in Aa(B) if the electron- 
electron interaction is not negligible in the diffusion chan- 
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29j . The corresponding correction to the conduc- 
tivity is Aai(B) — ^T^rg2(h) with h = Ez/ksT, where 
F a is a parameter reflecting the strength of dynamically 
screened Coulomb interaction. 

Since the electron g-factor of bulk Bi2Se3 is quite 
large [3(| , one would expect sizable Zeeman corrections to 
Aa(B). This appears to be in contradiction with the data 
recorded in tilted magnetic fields and plotted in Fig. 2(a). 
The low field magnetoconductivity exhibits very little an- 
gular dependence for tilt angles less than 80°. Consid- 
ering that E z nearly doubles (triples) for 6 =60° (70°) 
with respect to the zero-tilt case, we conclude that the 
influence of the Zeeman energy can be neglected in case 



of zero- or small tilts. In the non-interacting regime, 
this can be understood as a consequence of strong SOC, 
and hence small 7 for the fields of interest. Also in the 
regime where e-e interactions are important, the strong 
SOC suppresses the Zeeman contribution. The Zeeman 
term was derived under the assumption of weak SOC [28| . 
Theories |28l . I29I . [3lj and experiments [32j on other ma- 
terials have clearly shown that strong SOC can diminish 
and even entirely suppress the Zeeman-split term in the 
diffusion channel. 

The effects of strong SOC are further manifested in 
the temperature dependence of ct xx displayed in Fig. 2(b). 
The slope of the Aa(B)-\nT plot, defined as k = 
(Trh/e 2 )da xx {B, T)/d(\nT), is nearly constant for B=0.2- 
5 T. Both weak antilocalization and e-e interaction can 
cause the InT dependence [3 S|- The weak antilo- 
calization effect however only produces a pronounced 
T-dependence to cr xx at zero or low magnetic fields. 
The nearly constant slope at B > 0.2 T can be at- 
tributed to the strong SOC, which suppresses the triplet 
terms 0, |3l[ . They would otherwise produce In T cor- 
rectionsproportional to F a 28|, |3l| . Hence, K = 1 is ex- 
pected [5| for sufficiently large B. The observed n ~ 0.8 
is slightly smaller. This deviation may originate from 
other sources such as the corrections in the Cooperon 
channel (33|. Nevertheless, the nearly constant n for 
B = 0.2-5 T indicates that the e-e interaction does not 
induce significant corrections to Aa(B) in lower fields, 
where the weak antilocalization is pronounced. As to 
the zero-field conductivity, the combined effects of e-e 
interactions and weak antilocalization lead to a InT de- 
pendence with k = 1 — I/2=f/2, which is qual itatively 
in agreement with our data and others (Tol TTll| . Taken 
together with the tilted-field data, we can conclude that 
the single-particle HLN equation in the strong SOC limit 
(Eq. (1)) can provide a good description of the low field 
magnetoconductivity for a conducting 2D channel with 
strong SOC or one surface of a 3D TI. 

Now we are in a position to use the low field magneto- 
transport as a tool to analyze the influence of a nega- 
tive gate voltage. Fig. 3(a) shows Hall data from one 
of the samples with large density tunability. It is a 
10 nm thick undoped Bi2Se3 film. The electron density 
at Vg—0, estimated from the low field Hall resistance, 
is about 2.7 x 10 13 cm -2 . The Hall resistance, i? xy (B), 
increases as Vg decreases. It reaches a maximum at 
Vg = — I25V, which would correspond to an electron 
density of n e w 0.3 x 10 13 cm -2 . Further decrease in Vg 
leads to smaller R xy (B) and even reversal of its sign. For 
Vg < — 150 V, the Hall curves become strongly nonlinear. 
The high field Hall coefficient is plotted in Fig. 3(b) and 
depends non-monotonously on gate voltage. Also shown 
is the longitudinal resistivity at B — 0, denoted as p xx (0) 
throughout this paper. It also exhibits a non-monotonic 
dependence on Vg- This, together with the fact that the 
Hall coefficient Rh does not reach a minimum, points 
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FIG. 3. (color online) (a) Hall resistance curves for Vg = 
-125, -90,0, -170, -180 V (from top to bottom), (b) Gate- 
voltage dependence of p xx at B — and high field Hall coeffi- 
cient, defined as dR xy /dB and fitted from the data in B — 16- 
18 T. (c) Gate- voltage dependence of a and B$ obtained from 
fits to Eq. (1). Shown in the left and right insets are the band 
diagrams for large and small negative gate voltages, respec- 
tively. The top (bottom) surface is depicted on the left(right). 
(d) Vg dependence of B^i (hexagons) and B^ (triangles) ob- 
tained from fits to Eq. (2). B^i and B^ can be assigned to 
the bottom and the top surfaces, respectively. This 10 nm 
thick sample only has a back-gate. 



to the coexistence of electrons and holes for large nega- 
tive gate voltages. It is noteworthy that the maximum 
in p xx (0) appears at a Vg smaller than that of the Rh 
maximum. Therefore, the crossover from the pure elec- 
tron system to the electron-hole system must take place 
before the appearance of the /? xx (0) maximum. 

For the gate voltages smaller than that at the i? xy max- 
imum, the Fermi energy on the bottom and top surfaces 
are expected to lie below and above the Dirac point, 
respectively, even though the precise position of Ep is 
not known. As a consequence, the Fermi energy in the 
bulk (or at least part of the bulk) must be located in the 
band gap. The nearly one order of magnitude increase in 
Pxx(0) as Vg is lowered from to —150V is much larger 
than what has been reported for cleaved Bi2Se3 flakes 
cleaved on Si02/Si substrates The significantly en- 
hanced p xx (0) is an encouraging signature that much of 
the bulk conductivity can be suppressed. It can reach 
values as high as ~ h/e 2 0]. 

The magnetoconductivity also exhibits a strong gate- 
voltage dependence, especially for Vg < —50 V. Best fits 
to Eq. (1) yield the data plotted in Fig. 3(c). The most 
striking feature is that a is close to 1/2 for Vg > —70 V, 
and it increases to values close to 1 for Vg < —140 V. 
In the crossover region (—70 to —140 V), the Bi2Se3 thin 
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film undergoes a transition from a low density electron 
system to a separated electron-hole system. Hence, for 
large negative gate voltages, a fit of the magnetoconduc- 
tivity data to a two-component HLN equation is more 
appropriate: 



Acr(B) 



2 2 

—y 

2irh ^ 



(2) 



Here B^i and B^2 are dephasing fields for conducting 
components 1 and 2, respectively. As shown in Fig. 3(d), 
they have opposite dependencies on Vq. B^\(B^2) de- 
creases (increases) as Vg is lowered. They approach ap- 
proximately the same value for large negative gate volt- 
ages. The dephasing field is proportional to (Dt^ 1 cx 
(fj?T e T0) -1 , so B$ is expected to increase for the electron 
component on the top surface, while it should decrease 
for the holes at the bottom surface (interface) with de- 
creasing Vg- Therefore, the two curves with larger and 
smaller values of B^ in Fig. 3(d) could be assigned to the 
bottom and top surfaces, respectively. 

The observation that a increases toward 1 based on fits 
of the magnetoconductivity data to the single-component 
HLN equation (Eq. 1) implies that the top and bottom 
surfaces of the film make separate contributions to the 
conductivity [34] . Obtaining a values close to 1, how- 
ever, not only requires two decoupled conduction chan- 
nels, but also demands that both conduction channels 
have nearly identical dephasing fields. This is in general 
hard to achieve, in particular for samples where the gate 
tunability is not sufficient or the substrate surface is too 
rough (zj. Caution should also be taken to ensure that 
the transport is in the diffusive and weakly disordered 
(kpl 3> 1) regime for which the HLN equation is valid. 
For highly resistive samples, e.g. p xx ~ h/e 2 as shown in 
Ref. [7|, the condition UfI ^ 1 is no longer satisfied. 

In conclusion, it is possible to identify the surface 
transport of 3D topological insulators from magnetocon- 
ductivity measurements when transport takes place in 
the weakly disordered, diffusive regime. The use of a 
high-fc dielectric such as SrTiOs for back-gating enables 
us to tune the transport properties of both the top and 
the bottom surfaces. This device geometry is particu- 
larly useful for the future exploration of hybrid devices 
in which a topological insulator is interfaced with a su- 
perconductor HI, 1(| or a ferromagnet 15-17Q. 
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